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^ ■ The Schrodinger equation for stationary states is studied in a central potential V{r) pro- 

Qh! portional to r~^ in an arbitrary number of spatial dimensions. The presence of a single 

P3 ■ term in the potential makes it impossible to use previous algorithms, which only work for 

^ : 

^ . quasi-exactly-solvable problems. Nevertheless, the analysis of the stationary Schrodinger 

equation in the neighbourhood of the origin and of the point at infinity is found to provide 
/\ ' relevant information about the desired solutions for all values of the radial coordinate. The 

original eigenvalue equation is mapped into a differential equation with milder singulari- 
ties, and the role played by the particular case /9 = 4 is elucidated. In general, whenever 
the parameter (3 is even and larger than 4, a recursive algorithm for the evaluation of 
eigenfunctions is obtained. Eventually, in the particular case of two spatial dimensions, 
the exact form of the ground-state wave function is obtained for a potential containing a 
finite number of inverse powers of r, with the associated energy eigenvalue. 
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1. INTRODUCTION 

Many efforts have been produced in the hterature over several decades to study the sta- 
tionary Schrodinger equation in various dimensions with a central potential containing 
negative powers of the radial coordinate [1-19]. Some relevant examples are as follows. 

(i) When gaseous ions or electrons move through a gas whose molecules are not too large, 
then the two interact according to the law [3] 

V{r) = -^e'^ar-^, (1.1) 

where e is the ionic charge, a the molecular polarizability, and r the distance between the 
ion and the molecule. 

(ii) Interactions in one-electron atoms, muonic and hadronic and Rydberg atoms; photode- 
cay of excited states. 

(iii) Repulsive singular potentials, for which the stationary Schrodinger equation has non- 
Fuchsian singularities at zero and at infinity, may be sometimes expressed by a Laurent 
series [12-18] 



oo 



V{r)= J2 «-^" (1-2) 



-oo 



in the annulus r g]0, oo[. 

(iv) The Dirac equation for a spin-| particle interacting with scalar, electric and magnetic 
potentials leads to second-order equations for spinor wave functions involving an effective 
potential containing terms proportional to r~^, r~^ and r~^. Such equations can be used to 
prove the existence of magnetic resonances between massive and massless spin- 1 particles 
with magnetic moments [19]. 

Major progress has been possible when the potential contains a finite number of powers 
of r, e.g. 7--i^7--2 ,^-3 g^j^j^ j,-4^ Qj, ^2 ,^-4 ^^^^ ^-6 j-'j'j g^Q];^ g^ property is not a math- 
ematical accident, because in many cases one only succeeds in finding some eigenstates 



and eigenvalues for the so-called quasi-exactly-solvable problems [20]. What happens is 
that an algorithm can be found to evaluate some eigenfunctions provided that a number 
of coefficients in the potential remain non-vanishing and obey some restrictions, so that 
enough algebraic equations are obtained for all unknown parameters in the ansatz. 

The truly hard mathematical task, however, remains the one of solving the stationary 
Schrodinger equation with only one term in the potential, or for potentials admitting a 
Laurent series expansion for r g]0,oo[. The latter problem is studied by several authors 
(see [12-18] and references therein) and hence we focus on the former. When the stationary 
Schrodinger equation is studied for a central potential U{r) in R^, it reads [18] 
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(1.3) 



where /(/ + g — 2) results from the action of the Laplace-Beltrami operator on square- 
integrable functions on the {q — l)-sphere. Equation (1.3) can be re-expressed in a form 
which does not involve first derivatives, i.e. [21] 
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y{r) = 0, 



(1.4) 



where we have defined [21] 
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V{r) ^ ^U{r), 



A = / + ^(g-2), 



(g-l) 

y{r) = r 2 -i/;(r). 



(1.5) 

(1.6) 

(1.7) 
(1.8) 



having denoted by m the mass parameter, by E the eigenvalues, by / the angular momen- 
tum quantum number and by ip{r) the radial part of the wave function. The form (1.4) 
of the stationary Schrodinger equation is the most suitable for applications to scattering 
problems, because the operator in square brackets is an even function of A [12, 21]. 



For a given form of inverse-power potential, Sec. 2 studies the limiting form of Eq. 
(1.4) as r ^ and as r ^ oo to obtain a convenient factorization of its solution for all 
values of r. This involves an unknown function F, and an algorithm for the evaluation of 
F is developed in Sec. 3, if the inverse power of r in the potential V{r) is even. A class of 
ground-state wave functions are evaluated in Sec. 4, and concluding remarks are presented 
in Sec. 5. 

2. PROPERTIES OF EIGENFUNCTIONS 

For the reasons described in the introduction, we are interested in a potential V having 
the form (hereafter (3 > 2) 

V{r) = ar-f^, (2.1) 

where a and (3 are some dimensionful and dimensionless parameters, respectively. The 
first step in the attempt of solving Eq. (1.4) is the search of regular solutions in the 
neighbourhood of the origin. We therefore look for solutions y{r) taking the limiting form 

y{r) = r^e"^'^"' as r ^ 0, (2.2) 

where p, 7 and 6 are some parameters to be determined by consistency conditions, having 
taken S to be positive, as is suggested by known results for some values of (3, e.g. when 
(3 = 4. Indeed, on inserting the ansatz (2.2) into Eq. (1.4) one finds when r ^ the 
equation (for a fixed value of k, which is therefore negligible with respect to inverse powers 
of r) 



p(p-l)-(A2-i) ^ jd{2p-S~l) , 725 



2x2 
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^2 j,5-\-2 ^25+2 j,p ' V ■ / 

Since r ^ and 5 > 0, this limiting form of the equation reduces to 

ll!2lzAz21^J^-^o. (2.4) 

j,0+2 j,2o+2 j,p ^ ' 

Of course, the term proportional to r"^'^"^ dominates over the term proportional to r"^"^ 
as r — > 0, and hence Eq. (2.4) can only be satisfied if 

,2x2 



7^(5^ = a, (2.5) 



25 + 2 = (3, (2.6a) 

which imphes that the potential is repulsive (since a > 0), and 

5=^-1, (2.66) 

7=^^. (2.7) 

^ (/3-2) ^ ^^ 

So far, p remains undetermined. Note, however, that the coefficient {2p — 5 — 1) in Eq. 
(2.4) can be taken to vanish. This assumption leads to 

5 = 2p-l, (2.8) 

and by virtue of (2.6b) and (2.8) one finds /3 = Ap. In particular, if p = 1, one recovers the 
familiar limiting form of the solution when the repulsive potential is proportional to r~^: 

y(r) = re"^^"' as r ^ 0. (2.9) 

In the neighbourhood of the point at infinity, solutions of Eq. (1.4) behave as (hereafter 

e = ±l) 

y(r) ~ e^^'^^ asr^oo, (2.10) 

since only positive values of E are compatible with a repulsive potential. The ansatz for 
finding solutions of Eq. (1.4) for all values of r when the potential (2.1) is considered reads 
therefore 

y{r) = e"^'^"' e^^^^F(r), (2.11) 

where the function F interpolates in between the asymptotic regimes described by (2.2) 
and (2.10). By virtue of (1.4), (2.1) and (2.11) the function F should solve the differential 
equation 



J 2 +p(^):r + ^('') 

ar"^ ar 



F{r) = 0, (2.12) 



where 

p(r) = 2-f5r-^-'^ + 2ie^fK, (2.13) 



q{r) = 72,52^-2^-2 - ^5(5 + l)r-^-2 + 2ie-f6^/Kr-^-^ 



_«^-/3_(^l_i). (2.14) 



3. ALGORITHM FOR THE EVALUATION OF F{r) 

If we now assume that (2.5) and (2.6a) remain valid, the formulae (2.13) and (2.14) reduce 
to 

p{r) = 2v^r" 2 + 2ie^/K, (3.1) 



P r- -l-i .. , -I (^^ - ^' 

— \/ar 2 "^ +2ie\/aKr 2 _ 

2 r^ 



Q(r) = -^v^r-2-i+2i£v^^^~^ - ,^2^ " (3-2) 



Since we are taking /? > 2, such formulae lead to non-Fuchsian singularities in Eq. (2.12), 
and hence we look for F(r) in the form 

00 
F{r) = r'^ Yl ^'^' = ^""^(O- (3-3) 



s= — 00 



It is crucial to allow for negative powers of r in the solution, which are associated to the 
non-Fuchsian nature of the singular point at r = 0. By virtue of (2.12) and (3.3), the series 
a obeys the second-order equation 



(P f2u , \ d fuj(u;-l) cu , , , , 



a{r) = 0. (3.4) 



In Eq. (3.4), the term ^p{r) + q{r) has the coefficient ^/a (2a; — f j for r 2 1. We set it 
to zero to get rid of the dominant singularity at r = 0, which implies 

iv^^. (3.5) 

As a consistency check we point out that, if /3 = 4, one recovers the well known value 
uj = 1 (see (2.9)). For odd values of /3, uj is therefore a polydromy parameter for the wave 
function, as is familiar in singular potential scattering [13,17,18]. 
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If Eq. (3.5) is taken to hold, and if /3 is even, so that | is an integer fe, Eqs. (3.3) and 
(3.4) lead to the following recurrence relation among coefficients of the series a{r): 



2^/a{s + b + l)as+b+i + 2ie\^aKas+b 

+ 



(, + 2)(. + l)4^,(£-l)-(.^-i 



O-s+2 



+ 



2ie^{s + 1) + ^e(3^ 



tts+i = 0, (3.6) 



for all integer values of s from — oo through +oo. Such a formula is even more complicated 
than the recurrence relation for coefficients of Mathieu functions [22]. Nevertheless, it 
provides a well defined rule for the evaluation of F{r) for given values of «, when /3 is even. 
If /3 is odd, we no longer get power series, and hence we are unable to obtain a recursive 
algorithm. 

4. A CLASS OF GROUND-STATE WAVE FUNCTIONS 

It may be of some interest to conclude our paper by considering a different application of 
Eq. (1.4), i.e. the evaluation of ground-state wave functions for potentials V^(r) containing 
more than one negative power of r. In such a case an algorithm for bound states, rather 
than scattering states, can be found provided that the parameters in the potential obey a 
set of suitable restrictions, and we are now aiming to show how this can be obtained. For 
simplicity, we set q = 2 (bearing also in mind that two-dimensional models are of interest), 
so that A = / (see (1.7)), and we consider a potential in the form 

nr) = 4 + 4 + ^ + -- (4-1) 

ij^t: /yt'J ry^ /y 

The term -^ might indeed be combined with the angular momentum contribution -^^ — 3-^, 
but we prefer to keep them distinct to emphasize their different origin. Our ground-state 
ansatz for y{r) reads 

y{r) = expk{r), (4.2) 



where 



k{r) = — \- br + clog(r) a < 0, b < 0. 



(4.3) 



The negative signs of the coefficients a and b are necessary to ensure square integrabihty 
at zero and at infinity, respectively. The insertion of the ansatz (4.2) into Eq. (1.4) leads 
therefore to the equation 



y"{r)- k"{r) + {k'{r)y y{r) = 



(4.4) 



where the prime denotes the derivative with respect to the variable r. We now express 
y"{r) from Eq. (1.4) arriving at the equation 



-E + Vir)-^ = k"ir) + ik'ir)f, 



(4.5) 



because / = in the ground state. The above equations lead therefore to an algebraic 
equation where we equate coefficients of r^, for all p = —4, —3, —2, —1, 0. Hence we find 



a^ = A, 



(4.6) 



This system is solved by 



2a(l-c) =B, 



c{c-l)-2ab = C- - 



2bc = L>, 



r = -E. 



A, 



B B 

2a 2^/A 



6 = ^ 

2c 



D 



2 1 + 



B 



(4.7) 
(4.8) 

(4.9) 
(4.10) 

(4.11) 
(4.12) 

(4.13) 



B=-(,2 = , °' ,2 - (4.14) 

Note that, while (4.11) is consistent with a negative value of a as specified in (4.3), negative 
values of b are only obtained if D < and c > or if D > and c < 0. Moreover, after 
defining 

Eq. (4.8) can be expressed in the convenient form 

^ 1 . ^ D^/A 

C=_+^(l + ^) + ___. (4.16) 

For the evaluation of excited states, one has to write the ansatz for y{r) in the form of a 
product, and the resulting analysis is much harder. For this purpose, a separate paper is 
in order, which goes beyond the aims of the present work. 

5. CONCLUDING REMARKS 

The first original result of our paper is given by the formulae (2.11)-(2.14) and (3.1)-(3.6) 
for the solutions of the stationary Schrodinger equation when a central potential of the 
form (2.1) is considered in q spatial dimensions, if (3 is even and larger than 4. If /? = 4, we 
recover instead the standard result, according to which the regular eigenfunction behaves 
as re~^^ in the neighbourhood of the origin. Although much work had been done in 
the literature on similar problems [23], an investigation as the one we have proposed was 
missing to our knowledge. 

When the parameter /J is odd, or for non-integer values of /5 larger than 2, the equations 
of Sec. 2 remain valid, but it is no longer possible to develop a recursive algorithm as we 
have done in Sec. 3. At the mathematical level, the issue of self-adjoint extensions of our 
Schrodinger operators deserves careful consideration as well (cf. [24]). 

Last, in Sec. 4, the complete form of the ground-state wave function has been obtained 
in two dimensions when the potential takes the form (4.1). The energy eigenvalue is then 
given by Eq. (4.14), provided that the conditions (4.11)-(4.13) and (4.16) hold. 
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